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SY1T0PSIS 


The present work is an investigation into the dynamic 
response of a lathe bed. An experimental set-up was designed 
and fabricated for the purpose of studying the natural frequen- 
cies, mode shapes and the response locus of the lathe bed. 

lumped parameter technique is applied for analysing 
the bed structure. The lathe bed is modelled as a space frame 
with masses Imped at suitable stations and connected by mass- 
less spring elements. The analysis of bed is made (i) with 
port/al and (ii) without port/al. 

A computer programme is developed for analysis of the 
bed structure to get the natural frequencies harmonic response 
locus and the mode sh^ies. 

The dynamic response of the lathe was measured by 
exciting xhe structure at the centre by a magnetic shaker. P. 
force dynamometer was placed between the shaker and the struc- 
ture and the force signal from the dynamometer was displayed 
on an oscilloscope for measurement. Piezoelectric accelero- 
meters were used as vibration pickups and the signal from the 
accelerometers were displayed on another oscilloscope for 
measurement. 

The computed and experimental values were compared 
and a close agreement between the two was observed. 
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CHAPTER I 


INTRODUCTION AND REVIEW OF PREVIOUS WORE 

1.1 Introduction : 

The demand for better performance of machine tool re- 
quires consistent improvements of their structural design. The 
design requirements are all the more stringent in the case of 
programme controlled machines, where more exacting dynamic per- 
formances are required. 

Some attempts (Ref . 1 ) have been made for predicting 
the dynamic behaviour of machine tool structures using perspex 
model. However very little work has been done in the applica- 
tion of theoretical techniques to the actual machine tool struc- 
tures (Ref. 2, 3, 4). Very little use has been made of computers 
for predicting the dynamic behaviour of machine tool for select- 
ing the optimum structural details at the design stage. 

Vibrations produced due to variation of cutting force 
pose a serious problem in most machine tools specially when large 
volume of materials is removed at high cutting speeds. These 
vibrations result in poor finish of work as well as rapid tool 
failure. In general, machines resist vibrations better if na- 
tural frequencies of vibrations are high and has high damping 
capacity. The danping of a structure depends on (a) the 
materials from which it is built (b) the method of joining va- 
rious elements of structure and (c) machine foundations. 
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The vibrations encountered during operation of machine 
tools can be classified into three categories, viz. free, forced 
and self induced vibrations. The free vibrations are caused by 
mechanical impulse transmitted through the foundation. In fact 
the free vibrations die sway fairly quickly if the damping of the 
structure is high. Suitable floor mountings can be used to redu- 
ce the magnitude of vibrations transmitted through the foundations. 

Forced vibrations result from periodic forces either 
due to unbalanced rotating masses or caused by discontinuous cutt- 
ing processes, as is the case with milling operation. These 
vibrations grow to alarming levels when frequency of periodic 
forces is close to the natural frequency of the structure . These 
vibrations are dangerous even when first and second harmonics are 
close to the natural frequency of the structure. Forced vibra- 
tions are controlled by using balancing masses and by working at 
suitable spindle speeds. 

The self induced vibrations are due to the vibratio ns uTxV'ak- 0 *'' 4 
of the cutting forces caused by the cyclic -variation of the 
cutting parameter. These vibrations are called chatter vibra- 
tions. It is usually difficult to eliminate chatter vibrations 
in a given structure except by reducing the cutting speed. 

The self induced vibrations are mainly caused by 
dynamic instability of a machine tool. The .instability comes 
into picture because of the interaction between the dynamics of 
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cutting process and dynamic behaviour of machine tool 

structure. 

1.2 Review of the Previous work : 

Marten (Ref. 5) conducted experiments on actual structure 
and reduced scale model and concluded that the two did not give a 
close agreement of results. However he was able to spot the ,r eak 
points of the structure by experiments on the reduced scale models. 

Taylor and ^obias (Ref. 2) used Imped parameter method of 
analysis for _ analysing a milling machine. Further the computed values 
were checked by conducting the experiments on the actual machine and 
a fairly good agreement of the results was seen. 

Taylor (Ref. 3) analysed a planing machine on digital 
computer using the lumped parameter technique. Computed values showed 
a fairly good agreement with the experimental results. 

Gupta (Ref. 4) studied the dynamic response of a lathe bed 
using lumped parameter technique. He also carried out experiments to 
verify hie theor*tic*l analysis. However he did not get sufficient 
agreement between the computed and experimental results, lack of 
agreement was attributed to certain difficulties in experimentation 
and inaccuracies in modelling the structures. 
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The present investigation is a follow-up of Gupta’s work. 
The lathe bed analysed by Gupta (Ref. 4) is remodelled by evitably 
altering the end conditions and a computer programme is developed 
to calculate natural frequeneiesharmonic response locus, and mode 
shapes of the structure. Experiments are carried out to verify the 
theoretical analysis. 



CHAPTER II 


THEORETICAL ANALYSIS 

2.1 Introduction : 

In this chapter an outline of theoretical analysis 
is given; which is used for the analysis of the bed of a 
lathe. The various stages in analysis are : 

i) constructing e. geometrical configuration of 
lathe bed; 

ii) propose^ ft model with portals and ?l model with- 
out portals; 

iii) the coordinate systems ; 

iv) the mathematical simulation of structure; 

v) constructing the mass matrix and stiffness matrix 
of various elements; 

vi) assembly of stiffness and mass matrices of the 
whole structure; 

vii) the equations of motions and response of structure. 

2.2 Geometrical Configuration of the lathe Bed : 

The lathe bed is shown in Figure 1 and its photograph 
in Figure IS. The geometrical details of the bed were taken 
from Ref. 4. Eecause of non availability of detailed drawings 
the bed details were measured by a tape and recorded by Gupta. 
These dimensions were again checked and found to be correct. 
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2.3 Proposal f Model : 

Two models fox the lathe bed m&e proposed. Firstly, 
the bed with port^als removed and secondly the bed with port/als. 
Two models proposed such that the effect of po/trals on 

the bed behaviour could be observed; both theoretically and 
experimentally. The proposed models are shorn in Figures 1A 
and IB. 

In the model when port^sls are removed, the bed is 
bolted to the foundation by four bolts providing rigidly fixed 
end condition and in the model the beams are thus taken to be 
fixed ends. 

In the model when portjfels are there, po/tral is 
modelled by four vertical elements at each end and the ends are 
rigidly fixed to the foundation. Since we are using lumped 
parameter technique, the mosses have been assumed to be lumped 
at suitable stations. The stations selected are all the joints 
of rib and the beam and the fixed ends. By having stations as 
the joints it was possible to have the mass of ribs transferred 
on to the beams. The stations are connected by massless spring 
elements. 

2.4 The Coordinate Systems : 

. r, 

Two coordinate systems - local coordinate system of 
individual element and the reference coordinate for the 
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structure, are used in analysis. The local coordinate system 
has its origin at one end of the element and the length of 
element along its x axis. The reference coordinate system 
has its origin at the lower left corner joint of structure. 

The mass and stiffness matrix for each element is 
first written in local coordinate system and then transformed 
to the reference coordinate system for assembly. 

2.5 The Mathematical Simulation of Structure : 

The model with portfal, shown in Figure 1A has 22 
station points y and has seven types of elements. The station 
points can be numbered in an arbitrary falshon, but once num- 
bered the same numbering be referred to every time. The sta- 
tions are numbered as shown. This pattern is choosen since 
this gives the assembled mass matrix as diagonal one. The 
stations points 1, 2, 3> 4 and 19, 20, 21, 22 are completely 
restrained. The structure is modelled as a space frame, thus 
each end of an element has 3 rotational end 3 translational 
degrees of freedom. Thus the total number of degrees of free- 
dom for structure are 132, and of which 48 are restrained. 

The model without port*al, shown in Figure IB has 
total of 14 station points, from which 1, 2 and 13, 14 a 3 ?e 
completely restrained. The structure has 84 degrees of free- 
dom out of which 24 are restrained. 
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2-6 Construction of Mass and Stiffness Matrices : 

Por constructing the mass and stiffness matrices and then 

the equations of motions from them, a fev; assumptions are uade. It 

is 

is assumed that the bed /sade up of a few basic beam elements and 
elements are rigidly connected to each other at joints. It is 
assumed that the force displacement relations are unique, and are 
not (affected by path of loading the structure or motion causing the 
deflected shape. This means that material is ideally elastic over 
the range of stress and strain., involved in the loading; end the 
friction forces in structure and drag forces induced due to move- 
ment of structure in ambient fluid, air, are absent. Also the 
stress strain relations are considered to be the linear ones i.e. 
force and displacements are also linearly related. This allows us 
to make use- of the superposition principle. This assumption of 
linearity is highly idealistic and the force displacement relation- 
ship also depends upon geometry and constraints of structure. The 
linearity in force deflection relations cannot be postulated in 
ideally elastic materials, since the distribution of internal forces 
and hence strains change due to change in geometry as load and 
deflection are increased. It is also possible that during loading 
the system constraint may change thus introducing another source 
of non-linearity. 

However, we are working under sufficiently small dis- 
placements such that the force deflection relations are quite 
close to linearity. 
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The element stiffness matrix is made up cf elements which 
gives the force required in a given degree of freedom to produce a 
unit deflection in the other given degree of freedom of element. 
There are veil known standard techniques for obtaining the elements 
of stiffness matrix (Ref. 5) vis. - by inversion of force displace- 
ment relations, unit displcc. lent theorem, Castiglanc’s theorem. 

A brief outline of the procedure to obtain the stiffness matrix is 
given in appendix II. The sections! properties cf bean required 
in computation cf stiff matrix element are calculated using the 
procedure given in Ref. (6 ) and is given in appendix I. 

The element mass matrix is made up of mass and inertia 
terms. The sectional properties of elements are made use of for 
the calculation of mass. The density for the material is taken 
from Ref. ( 4 ), Reference ( 5 ) gives a standard method of obtain- 
ing the mass matrix. A brief account of it and of the form of 
mass matrix is given in appendix III. 


2.7 Assembly of Stiffness and Mass Matrices : 

The element stiffness and mass matrices obtained are in 
the elemental coordinate. These are first transformed into the 
, reference coordinate by a suitable transformation defined by 


! T ii T \i*< T ~ - i-' r i -1 

•; k 4 = I A j £k jf^l 


( 1 ) 


where f" k 1 iis the element stiffness matrix in the reference coor- 
dinate system, j k^ being inlocal coordinate system an&jj^Jis the 
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transformation matrix obtained by resolving datum displacements 
in the local coordinates directions. Appendix (TV) shows that 
^he [ %J 1- s dsde up of direction cosines of the angles between 
local coordinate end reference coordinate system. The fora of 
is 

r \. 
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( 3 ) 


and 1, m, n are the direction cosines of the local coordinates ox, 
oy, oz in the reference coordinates ox, oy, oz-. 
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Tiie combined stiffness matrix is reached at by expressing 
the elemental displacements in terms of displacements of structure 


I d * [d['’] ^ 

in which £ ujis the elemental displacement in reference coordinate 
system and £. T lji is structures displacement in reference coordinates. 
£ A J is square matrix having elements as 0 and 1 and thus relating 
elemental displacements tc structures displacements. If is 

the generalised load vector* i.e • 


O] 

where P. is lead in direction cf diaplaeement u.., 

[ = f ¥ 1 k 2 k 1 k*J 


(5) 


and writing 


where£k^|is a matrix having) k 1 | submatrices at its diagonal. 
Now give a vertual displacement f S u i tc the system and relating 
the force and displacements we have (Hef .5 ) 


[P] - [AfOJtAjtSJ (6) 

Also by force displacement relation 

[P] = (7) 

Thus [I] = £AJ* r : ]p ] 

Equation 8 gives the combined stiffness matrix. Equation 7 gives 
a static equilibrium condition when the load matrix Dd includes 


the reactions as well. Since all these fo.rces are not independent 
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because of the overall equilibrium of structure, the matrix 
is thus singular . The six dependent equations corresponding tc 
six degrees of freedom of rigid body motion are to be eliminated. 
This is done by making the six displacement on certain selected 
point equal to zero and eliminating corresponding rows and columns 
from the£x J . This gives the following relation 

[ p ] r = [ r\l (9) 

r indicates that these are reduced matrices. 

Equation 8 can also be written as 

JW- 


This summation is easier to perform rather then computing the 
matrix multiplication. The boundary conditions can easily be 
incorporated here in the assembly. The same is achieved on com- 
puter by using a chart called Assembly Control Chart, which is 
no -thing but a simple table giving details of geometry of structure. 


i K 
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I f 
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i=i 




(10) 
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The assembly control chart in fact adds up all the elements of 

elemental stiffness matrix coming cn to the particular location 

* 

in the complete matrix. 

Proceeding on the some lines as above the mass matrix 
is also assembled. Pirst it is converted from local coordinate 
system to reference coordinate system by CXI as under 

C = 3-W*C»Jt>J 

and by a similar procedure as in case of stiffness matrix the 
mass matrix of complete structure is 

[IJ- (12) 

Again this assembly can be done by summation and the matrices 
be reduced as in case of stiffness matrix. 

2.8 Equations of Motions 

The equations of motion for the structure, when a perio- 
dic force cf amplitude P^ and frequency CO is applied, are 

+ C c 3 W = [p 0 3sin.wt. (13) 

where Jis the damping matrix £ Jis the displacement vector. 

is 

The value of dsnping constant/taken from Pr.f . & . This value 
was checked by striking the bed end storing the signal picked by 
accelerometer on a storage scope. Other investigators have said 
that the damping is proport ionr.l to mass matrix or stiffness 
matrix or a linear combination of the two (Ref. 2, 3, 7). 
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To solve equation fl3) for response end natural 
frequencies, we proceed as under {Ref. 5). Define a trans- 
formation matrix £ N J such that 

C Nj = CmJ " 1/2 (14) 

and W=[sJ[Y] (15) 

Then the equation of motion can be written as 

+MWltl 1 = QP 0 lSin-:t (16) 

Premultiply (16) by I 



\ SQNjfY 


.T 


[IT ; [ G j,‘ If ' Y | +[ F f i K ; 1 N { j Y/ 


il^PolSin. 


because ( 14 ) 


( 17 ) 


}_&M Mjflj = jl] 

denoting [ !TJ T [Gl[N] by [‘A] , fllf \ Yjll'by £b] and [ IT j T [ P j 
by { P^ j f then we have equation ( 17 ) as 


Wro + [^CYj+ c?JC Y J =C p il Sin ^ it (is) 

low we define another matrix transformation, such that 

W = M Lzj ( 19 ) 

and £ 3 ] is a matrix simultaneously dirgcnslising I A j and i B ] 


and 


L E f C B 3= tXJ 


( 20 ) 
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Equation (IS ) then changes tc 


3 3 + [ A J' H) L Z> r „B 1 \ H j l Z j = ■ P^j 


Sinn t . 


Prenultiplying by ! Rj we have 


« Hi [.I I ;R .1 [Z J + ’..E . L, iSi‘ Zj + ( Rj V BjUOfZ . =iH) jP^iSi 


Sinc^t 


Since\_ Aland fBj axe simultaneously being diagonalized , hence 

Lntz'3 +£:eJ Cz] +MW ; [t ] sm fe 

where 

r ■ ® ' j - L*f [ A] [.B J = [' 2 4 !« ni J 

I" *-j =Ih] t Ib;i|e ] -f» nl N J 

sndiP.1 = iR; T t P,} 

Since and are diagonal matrices, equation (22) repre- 

sents a set of uncoupled equations which can be solved for Z. 

One of the equation of the set is 
** - 2 

Z. + 2 <^ . Z . + ..a . Z. = P. sin c. t 

3 nj 3 nj 2 3 

and has a solution as 


Z. = 
3 


P . Sin co t 
3 

2 \2 , / , <2 2 


r / 2 2 \2 y.j 2 2 2v 1# 

L ( "nj 5 + ( 4 ^i“ni J 


at a lagging phase of 0 given by 


0 = Tan ( 


! Vi <^n.i 


(24) 
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Once these values of Z's are obtained these can be substituted 
bach to get the response in X, by using equation (19) and (14-). 


2.9 Frequency Computations : 


The natural frequencies are calculated for the system 
by obtaining the Eigen values of the matrix \ 3 ] . The Eigen 
values were obtained by Jacobi's method. In this method the 
matrix [.B \ is successively transformed to a diagonal matrix } 
using transformation matrices l RVs such that we have 


[ \ Vi v 2 -- % Bj /. 1 1 * ** ** - r » 

and then the Eigen values are the diagonal terms of matrix jpl>.j • 
The Eigen vectors corresponding to these Eigen values are, then 
the columns of the matrix *i T J given by 

1*3 M.*, 1 h 1 *,* — V 1 

These rotational matrices T B ’ are formulated on the 
basis of eleminating largest off diagonal term (Ref. 8). Since 
the Jacobi's method starts.with .the dominant Eigen value and we 
are interested in lowest Eigen value, the matrix (Bj is read 
in with a negative sign giving us the lowest Eigen value. 
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EZPEREiEMT.'i SET - T JP /US TEBTSMEmi-TTQS 


3.1 Introduction : 

To verify theoretical analysis the dynamic response of 
the structure was experiment ally determined taking lathe bed 

(i) with portral (ii) and without pcrtaral. The experimental set- 

/'* 

up is shown in Figure 2/ and the block diagram of instruments is 
shown in Figure 4. The lathe bed was excited at the centre/by 
an electromagnetic shaker and the force and displacement were 
measured using a force dynamometer and accelerometer respectively. 

The phase difference between force and the displacement was mea- 
sured with respect to a standard signal. 

3.2 Set - Up and Instrumentation : 

The electromagnetic shaker was held over a mild steel 
plate which was rigidly fixed to the foundation by six foundation 
bolts. The shaker generates a harmonic force signal. The shaker 
is connected to the dynamometer by a left hand - right hand bolti F *B , ^J 
The electrical power to shaker is given through a suitable pow/er 
amplifier. The power -amplifier operates on .2 mains . 

The force dynamometer (Fig. 3) was connected to shaker 
from one side and to the lathe bed from the other by suitable 
bolts. The dynamometer was designed [for a. load level of 50 lbs. 
amplitude of a periodic force, with a natural frequency cf 
1800 c.p.s. It was expected that in the range of 50 to 600 c.p.s 
the dynamometer would give a flat response. The dynamometer is 
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capable of measuring two forces s imult aneo us ly , mutually perpen- 
dicular as shown in Figure 3 ^22- Strain gauges we re pasted at 
suitable positions cn the dynamometer to minimise the ■ cross 
response. The dynamometer was calibrated in the two directions 
end the cross response was found to be negligibly small . The 
output from dynamometer was fed to an amplifier sc that the signal 
level could be raised tc a value suitable for the phase meter. 

The output frem dynamometer was also displayed cn an oscilloscope 
for measuring the level cf forcing function from shaker. 

The lathe bed was excited at the central rib by gripping 
central rib in a U belt. The TJ bolt was connected tc the bolt 
from shaker through a suitable plate . 

The Endevco accelerometers were used as vibration piek- 

, ( f i? z 0 

ups. Suitable I luminium blocks/ were mounted on the lathe bed 
and the accelerometers were mounted on these Aluminium blocks. 

The output from tin. accelerometers were amplified on a suitable 
amplifier and then displayed cn an oscilloscope and was also 
fed to the input terminals cf a phase meter. 

A phase meter was used tc- measure the phase difference 
between the forcing signal from dynamometer and the accelerometer 
output. Since the two signals were not upto the required level 
for the phase meter, they were compared with another standard 
signal of the same frequency and the phase difference was then 


computed 
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The frequency of excitation was controlled by a frequency generator. 

The output of frequency generator was fed to the power amplifier 
supplying power to magnetic shaker. 

3.5 Experimental Procedure : 

The experiment was first performed without portal with the 
bed bolted on the foundations. The following procedure was used: 

(1) The dynamometer was caliberated using standard weights. 
The calibration curve is found to be linear and to shown 
in figure 5. 

(2) The function generator was set to give sinusoidal wave 
form at a particular frequency. The output from function 
generator and power amplifier were checked on an oscillo- 
scope to ensure that the signal supplied to the shaker 

in sinusoidal. 

( 4 ) The power to the shaker was supplied by the power 
amplifier and amplitude of the forcing signal was 
maintained at 7.5 kgs. by adjusting the input current 
to the shaker. 

( 5 ) The signals from the accelerometers and the dynamometer 
were measured by displaying them on the oscilloscopes. 
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(6 ) The output from the dynamometer was fed to the phase 
meter through an amplifier and its inphase and quadrature 
components with respect to a reference signal from the 
function generator were measured. 

( 7 ) The inphase and quadrature components of the signal from 
the accelerometer were measured in a similar manner as 
under (6) by replacing the dynamometer signal with the 
accelerometer signal. 

(8) The above procedure was repeated by increasing the 
frequencies from 60 to 1000 cps in steps of 10 cps. 

Whenever a mode was detected, a number of readings in 
that range were were taken at suitable intervals. 

(9) The above experiments were repeated by fixing the portals 
between the bed and the foundations . 

The measured accelerations were converted to displacements by 
2 

dividing by . The phase between the forcing signal and the displacement 
was calculated by taking the difference of phase angle measured under 
(6) and ( 7 ) and by making appropriate corrections for the fact that the 

O 

accelerations signal lags behind the displacement by 180 . The phase 
shift in the amplifier used for amplifying the force signal was measured 

at various frequencies and is shown in figure 6. The average phase 

o 

shift of 6 was also accounted for in calculations. 
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(6) The output from the dynamometer was fed to the phase 
meter through an amplifier and its inphase and quadrature 
components with respect to a reference signal from the 
function generator were measured. 

( 7 ) The inphase and quadrature components of the signal from 
the accelerometer were measured in a similar manner as 
under (6) by replacing the dynamometer signal with the 
accelerometer signal. 

(8) The above procedure was repeated by increasing the 
frequencies from 60 to 1000 cps in steps of 10 cps. 

Tflhenever a mode was detected, a number of readings in 
that range were were taken at suitable intervals, 

(9) The above experiments were repeated by fixing the portals 
between the bed and the foundations. 

The measured accelerations were converted to displacements by 
2 

dividing by . The phase between the forcing signal and the displacement 
was calculated by talcing the difference of phase angle measured under 
(S ) and ( 7 ) and by making appropriate corrections for the fact that the 

O 

accelerations signal lags behind the displacement by 180 . The phase 
shift in the amplifier used for amplifying the force signal was measured 
at various frequencies and is shown in figure 6. The average phase 

O 

shift of 6 was also accounted for in calculations. 



CHAPTER 4 


HESUETS Alii) BISCUSSIC® 

4.1 Results : 

The frequency response curves, mode shapes and harmonic 
response locii have been plotted for computed and experimental values 
in Figures 7 to 17. 

The following constants were used in computation. (Ref. 4) 
Youngs Modulus E = 9.6 x 10^ p.s.i 

Poisson f s Ratio = 0.27 

Damping Coeff. = 0.0046 

For the bed with portal, frequency response curves at 
central rib have been plotted in Fig. 7, frequency response curves at 
first rib for experimental and computed values have been plotted in 
Fig. 8 and the direct harmonic response locus at central rib in 
vertical direction for experimental and computed values have been 
plotted in Fig. 9. For the bed without portal, frequency response 
curves at central rib have been plotted in Fig. 10, frequency res- 
ponse curves at first rib have been plotted in Fig. 11 and the direct 
harmonic response locii at central rib in vertical direction have 
been plotted in Fig. 12. The frequency response of foundation has 
been plotted in Fig. 13. Figures 14 and 15 show the experimental 
mode shapes without and with the portal, respectively. Figure 16 
and 17 show the computed mode shapes without and with the portal, 
respectively. . « 

4.2 Discussion : 

1. Bed with portal - The frequency response curves in 

Fig. 7 Show three i Bn! nr> emen+ r\e oVq mVio rvoolre o+ - A T ' ** 
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200 c.p.s and 290 c.p.s on computed curve with the major mode being 
at 290 c.p.s, and at 95 c.p.s., 130 c.p.s., 180 c.p.s. and 310 c.p.s, on 
the experimental curve with the major mode at ISO c.p.s. . The dis- 
placement at 310 c.p.s. was measured at a force level of 3.75 kg. 
amplitude, since it was not possible to excite the structure at 
7.5 kg. amplitude force level in the range of 300 to 400 c.p.s. .The 
measured displacement was doubled to correspond to the force level 
of 7.5 kg. amplitude. The errors in computed displacements at peak 
points are - 45 $ at 290 c.p.s., 8.45$ at 200 c.p.s. and 15*8$ at 
130 c.p.s.. The percentage errors in the theoretical natural frequen- 
cies are 11.1$ at 180 c.p .s.in verticel direction, 4.15$ at 230 c.p.s. 
in transverse direction and 3*53$ at 410 c.p.s. for torsional vibra- 
tions. The response locus curves (Fig. 8) show three distinct loops 
corresponding to three frequencies of peak displacements which is 
evident in frequency response curves (Figs. 7? 8, 10 and 11 ). 

2. Bed without portTal - Frequency response curves (Fig. 10) 
show three displacement peaks. The peaks occur at 335 c.p.s., 230 c.p.s. 
and 130 c.p.s. on the experimental curve and at 350 c.p.s., 225 c.p.s. 
and 160 c.p.s. on computed curve. The errors in displacement at peak 
points are 42$ at 350 c.p.s,, 14.4^ at 230 c.p.s. and 5.5$ at 130 c.p.s. 
The percentage errors in natural frequencies are 1*£ at 228 c.p.s. 
in vertical direction, 6$ at 275 c.p.s. in transverse direction and 
2 . 25 $ at 450 c.p.s. in torsional vibrations based on the theoretical 
values. The response locus curves show three distinct loops 
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corresponding to peak displacement as shewn in frequency response 
curve s also . 

Figures 7, 10 and 11 show that computed frequencies are 
higher than the experimental ones, ^hr lower natural frequencies 
obtained from experimental data can be explained on the basis of 
decreased rigidity of the supports at the ends as compared to the 
theoretically assumed end condition of complete restraint. The 
loss of rigidity brings down the natural frequencies and experi- 
mental values of displacements are larger. 

The study of mode shapes (Fig. 14 to 1?) reveals that 
the bed is weak: in vertical direction as the first mode is observed 
in this direction. The seflond mode corresponds to transverse direc- 
tion and the third one to torsional vibrations. The effect of por- 
t#al with bed (Fig. 14 to 17) is to lower the natural frequency due 
to increased flexibility. 

The reasons for discrepancies in displacements and frequen- 
cies observed experimentally are : 

i) Flexibility of the exciting system - the exciting 
system constituting the dynamometer, two bolts, the 
TJ bolt and the plate has its own natural frequency 
which affects the excitation and hence response of 


the structure . 
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li) flexibility of the supporting system - Due to 

the flexibility of the supports the boundary con- 
ditions considering fixed ends are not strictly 
true. The flexibility in the transverse direction 
changes by an appreciable amount due to the bolts 
at the supports. This lowers the natural frequency 
affecting largely the first mode. This flexibility 
also gives rise to rigid body motion, thus affect- 
ing the response. 

iii) ?eed back loop formed through the foundation also 
changes the dynamic characteristics of Hie bed. 



CHAPTE? 5 


OOFCI'USIONS /IIP SUGGEST ICES ?03 F T IRTEFF ' r C?2{ 

From the results of the present T .ork it is concluded that 
the lumped parameter technique when applied to the lathe bed can 
predict the natural frequencies within ICf? in the worst case, fur- 
thermore this technique can be used to compare two or more machines 
for their dynamic behaviour and it can also be used for selecting 
the best design out of the available sections at the design stage. 

For further work it is suggested that the structure be 
analysed along with mountings on the bed and response at the tool 
post, head stock and tail stock be measured. It is suggested that 
a greater number of vibration pick-ups be arranged and used such 
that the response at a, number of places be measured. Filters need 
to be used with these pick-ups such that law frequency signals and 
noise can be eliminated. It is suggested that to avoid the feed 
back loop the structure be excited from the top. Investigations 
can also be carried out by pre-loading the structure . 

It is also suggested that the structure be analysed with a 
greater number of station points. 1 further step would be to con- 
sider the distributed mass matrix of elements instead of lumped 
mass system and the damping matrix be estimated theoretically. 
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APEEEDK - 1 


Evaluation of Sectional Properties : 


the section of the element is approximated bp straight 
line sigments and circular arc segments* The procedure is to cal- 
culate the sectional properties under the upper contour and subtract 
the same of the lover ecntsrer of section, the fcrmulaes for t 
straight lines and circular arcs are given as under - 

(l ) Straight line 
Let a = (x 
b = 


2 X 1 } 


then 


A = * 

X 



1 


A is the area under line AB. The first moment of inertia and second 
moment of inertia about x axis and y axis and product moment of 
inertia are given as 


2 

2 


X 2 


2 j 

7 dx 




= \ s. (y.j (y^ + b ) + ) 


(a— 4) 


1 

2 


X 2 




xy dx = a (x^ (y^ + 'f ) + \ 2 tp + \ *>)) (a-5) 


I = 


1 

3 \ p dx = J a (v^ (y^ + -| b) + b^ (y 1 4- •£)) (a-6) 

X^ 


I = 


2 / 2 / b\ / a \ 

v dx = a (x 1 (y + - ; + a y 1 (x^ + ~ ) + 


|b a + |a )) 


(A-7) 
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"1*^2 i £ i 

I xy = 2 j to = - a (y 1 (x 1 + ^ a) + b ^ 

(y, + | ) + fab ( ?1 + fb)) (-4-8) 

(2) Oircular ire 

let © 1 = Sin 1 (y 1 - ^ ) (--.-9) 

6 2 = Sin" 1 (y 2 - | ) (A-IO) 

Thus for a convex arc 9^ ^ 9 p 
and 0^-9 77' 

and for a concave arc © ? and T 7 " <' © :< 2 to 

The area /, first moment of section S s and moment of inertia of 

section about various axes are given by 

/' X 2 2 

; y dx = ( 2 ( ) (Cos © 2 - Cos © 1 ) + 

X 1 

(Sin 2 © 2 - Sin 2 ©^ ) - (© 2 - ©^ )) (A-10) 

\ f 2 / <te = ((( | f * | ( cos S 2 - 0 ob e n ) 

x i 2 ' 

+ j ( -|) (Sin 2 © 2 - Sin 2 6^ ) 

"*12 ^ CoS 3 e 2 “ CoS 3 ""I ^ 6 2 “ ^1^ 

(A-1 1 ) 

X 2 

S = ] j xy dx = R 3 (( | ) ( | ) ( Cos © 2 - Cos ©.,) 

1 (Sin © 2 - Sin ©^ + | ( | ) 

(Cos 2 © 2 - Cos 2 ©^ j (| (Sin 2 ©g-Sin 2©^ 

+ ^ (Sin 3 © 2 - Sin 3 © 1 ) - (« 2 ~ © 1 )) 

(A-12) 
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I = 
x 


I = 

y 


i 

xy 


1 

7 y ax 

x i 


= E 4 ( 7 (~r + j ) (Cos Q 2 - Cos © 1 ) 

3E 2 

1 h 2 1 

+ 7 (-% + 3 ) ( Sin 2 g 2 ~ Sln 2 ®i ) 

H (Cos 3 © 2 - Cos 3 © 1 ) - y 


96 


1 h 1 

(Sin 4 © 2 - Sin 4 © 1 )- 2 ^ J? + 4 ^ 


( ®2 " 9 i^ 


(.4-1 3 ) 


i 2 . 

‘ y x dx 


1 

2 


H 4 ( | ( % + | ) (Cos 0 2 - Cos e 1 ) 

-i ( f ( (s “ V s±ne i ) + la ( i 

1 2 

(Cos 2 © 2 - Cos 2 0 ) + (Sin 2 © 9 

R 

Sin2 9^ +-|| (Cos 3 © 2 - Cos 3 0.,) 

1 2 h 

+ j — (Sin 2 6 0 - Sin 2 ) + 


1 7 ‘ 12 R 


(Cos 3 © 2 - Cos 3 0 1 ) +g§- (Sin 3 © 2 

- Sin 3 © 1 ) + (Sin 4 © 2 - Sin 4 6^) 

- ( — 2 + 7 ) (© 2 - © 1 )) (A- 14 ) 


' ^ ^ ( “S { T2 + f > (C ° S ®2 " C ° S 8 1 

“ 4E ^ Sin 9 2 “ Sin e i ) + Q ( “2 + 7 ) 

R 

(Cos 2 © 2 — Cos 2 © 1 ) + 1 &- 

(Sin 2 © 2 - Sin 2 ©^ - (Cos 3 © 2 
' C0S 3 ®l) + 11 (Si» 2 © 2 - Sin 3 ® 1 ) 

- h (Cos 4 ®2 - Cos 4 V * 

2R 

(A-1 5 ) 


(© 2 - ©.,)) 
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If the superscript 1 denotes the upper boundary and 
2 denotes the lower boundary of section thus for section we have 


A = 

- 


S = 

X 

s™ - 

X 

X 

S = 

y 

- 

y 

S. (2) 

J 

i 

X 

i (l > - 

X 

I< 2 ’ 

X 

I = 

y 

- 

X 

,(2) 

X 

c<J 

II 

- 

xy 

I< 2 > 

xy 


(A-16 ) 


Transformation to principal axes : 

The angle of rotation of coordinates is given by 


•, 2 I I' 

8 = 1 Tan -1 ( S I-r ) 


- 1 * 

t T 

where 1^ and 1^ are given by 

I ' = I Cos 2 9+1 Sin 26-1 Sin 26 
xx y xy 

I ' = I Sin 2 6+1 Cos 2 6+1 Sin 2 6 
y x y xy 


(A-17) 


(A-18) 

(A-19) 


Inertia Properties : 

The mass moment of inertia of the elements are obtained 


by relation 
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mxx 


2 2 

(x + y ) dx .dy .dz 


. , 2 - 2 
= .' { z dy.dz.dx + < ' y dy « dz • dx 


:• ( £. , _5_ 

> V 2 2 


5 ■ — ( V + 1 . ) 


where K = mass/unit volume = specific weight / g. 


nyy 


' 2 ^ 

j (z + x ) dx . dy . dz 

'.y t -z 

A 1 


I , ,3 - , ,3 

f ( I 4. — „„ \ = — , / T j- 

v z 24 ; 2 v y 


a r 
12 


and I 


mzz 


= r 


2 2 

(x + y ) dx . dy . dz 


> (- 


a r 

24 


’ P , -.3 


(A— 20) 


( A-2 1 ) 


(A-22) 
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Stiffness Matrix : 

The complete derivation of stiffness matrix .‘is given 
in reference 5 . 

The element k. . of stiffness matrix is defined as "the 
force required in the direction 3 to produce a unit deflection 
in the direction i." 




U*-J 


AHPEEDIX 3 


Mass Matrix : 

Che complete derivation of mass matrix is given in 
reference 5. The elemental mass matrix for the lumped parameter 
model is given by : 


Pi 




HDQC 


niyy 


m 


3 = 


mzz 


el 


m . 


el 


HDCX. 


myy 


mzz 



APEEIJDIX - 4 A 


Direction Cosine Matrix : 


'Xr- x -i 


Z-X-T, 

f 

\ | 



/ 

/ 

/ 

.TC 


A point P in space has coordinates (x, y, z) along ox, oy 
Por the new axes, ox, oy, oz, the coordinates of point are (x, y, z 



To relate 

(x, y, 

z) and ( x, y, z ), we have 

*1 

= hi x i 

+ X 12 

X 2 

+ 

1 13 X 3 

X 2 

= hih 

+ 1 22 

X 2 

+ 

X 33 X 3 

*3 

” hi 1 ! 

+ X 32 

X 2 

+ 

1 33 X 3 

where 1 

. . = Cos 

/ ^ 

(x . x . 
i 3 

) , 

the 

direction cosines. 


Expressing in matrix notation 


I x j = iAJfxJ 


where 


li 

■ ‘/ w ****’h 

IH* 

"x 1 

li 

[xl 

“ | X 1 


x i 

» 

; x 2 


-h] 


i yr 

1 

11 

’hi 

1 12 

1 13 


X 21 

1 22 

1 

23 


hi 

X 32 

X 33 j 


j oz. 
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Transformation of Elemental Matrix from local Coordinates to 
Reference Coordinates : 


The displacements in Ice al and reference coordinates are 

related by 


L U J = lAj 1. u J 

Introducing a virtual displacement, u in reference coordinates 
can be written in terms of u as 

\ bu J = U ] to u j 

If the loaning vector is P on local coordinates and P in 
the reference coordinates the Equating Virtual Work in reference 
and local coordinates we have 


or 


l-u] \p] = '(Wf [P] 

1 _ , T 

[Su] Ip] - jsa 1 t p j = L 0 J 

T T T 

[■bu j L P j - ^uj f;,] f PT = V 0 j 

_ T 

Lp)= ‘ ;:0 L PJ 


Since 

and 


[P J = ’..k 1 j u j 

1 P J = L k ' 3 t u j 


Hence i k x 3 L u i 

L k i ) f u j 
_ i 

or [k ! 


1 - _ 

' ( "i i _ 3L \ « Hi 

= 1 k j L u -> 

T 

= IaT ij 

= ur 
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DYNAMOMETER DE5IGIT : 


The two force dynamometer design is based on expressions 
given in reference (iC), The criteria fer design is the natural 
frequency of the dynamometer which should be 4 tc 5 times the fre- 
quency of forces to be measured, so that the recorded force is not 
influenced by any vibrating motion of the dynamometer, por the 
analysis purpose the dynamometer is reduced to a mass, m, supported 
by a spring. The natural frequency of such s system is given 


,_o = 

n 


_L I I 

2rr J m 


c .p.s 


where k is the stiffness of spring. 


In terms of weight of dynamometer w 

\ 

K 

‘"n “TV i 586 I °- p ' S 
K is in lbs ./in. and W in lbs. 

The stiffness of dynamometer in radial direction and tangential 
direction are 


K = 


E b t 

1.8 r 


and. 


L = 


E b t' 


3 .6 r 

where E is the modulus of elasticity of material of dynamometer, 


b is the width, t is the minimum thickness and r is the radius of ring 
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The strains encountered are given by 

,c ^ 10.9 F . r 

ir A = ± I 

a 2 

Ebt 




The above equations were used for getting the various dimensions 
of the dynamometer. 
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